arXiv:1505.03681vl [cs.CG] 14 May 2015 


A light metric spanner 


Lee-Ad Gottlieb 

Department of Computer Science and Mathematics 
Ariel University* 

May 15, 2015 


Abstract 

It has long been known that d-dimensional Euclidean point sets admit (l-l-e)-stretch spanners 
with lightness We = that is total edge weight at most We times the weight of the 

minimum spaning tree of the set [DHN93) . Whether or not a similar result holds for metric 
spaces with low doubling dimension has remained an important open problem, and has resisted 
numerous attempts at resolution. In this paper, we resolve the question in the affirmative, and 
show that doubling spaces admit (1-|-e)-stretch spanners with lightness We = (ddim /£) 0 (ddim)^ 
Important in its own right, our result also implies a much faster polynomial-time approxi¬ 
mation scheme for the traveling salesman problem in doubling metric spaces, improving upon 
the bound presented in |BGK12j . 
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1 Introduction 


Let G = {Vg,Eg) be a metric graph, where vertices Vg represent points of some metric set S, 
while the edge weights of Eg correspond to inter-point distances in S. A graph R = (V/j, Er) is a 
(1 -|- e)-stretch spanner of graph G = {Vg,Eg) if -R is a subgraph of G (specifically, Vr = Vg and 
Er C Eg), and also dR{u,v) < (1 -|- e)dG{u,v) for all u,v € G. Here, dG{u,v) and dR{u,v) denote 
the shortest path distance between u and v in the graphs G and R, respectively. 

Low-stretch spanners have been the subject of intensive and broad-range study over the past 
three decades. Research has focused on minimizing such properties as construction time, vertex 
degree, graph diameter, and total edge weight - as well as possible trade-offs between these quan¬ 
tities - in various settings such as planar graphs or Euclidean spaces [Vai9H ISal9H ISoa94( IAMS94( 
I(1K95[ rAWl IDN971IAMS991 ^^^+9"^ IDNS95[ ICTLNfM iROMlUl IDESORl lESHl IDHN931 \Mu] . 
Of particular interest is a remarkable result of the nineties, that d-dimensional Euclidean spaces 
admit (1 -|- e)-stretch spanners with lightness We = meaning that the total spanner weight 

is at most a factor We times the weight of the minimum spanning tree (MST) of the set [DHN93) . 

The work of Gao et al. [GGNOb] was the first to consider low-stretch spanners in metric spaces 
of low doubling dimension ~ a strictly more general setting that Euclidean space. This immediately 
spawned a long and fruitful line of work, showing that results comparable to those of Euclidean 
space can be obtained for these spaces as well, for construction time, degree and diameter. [CG061 
IGROSbt IGRn8a( ISmin9( IGKK131 IGLNS15( LSoll4| . However, the proof of lightness for Euclidean 
spanners relied heavily on the properties specific to that space (in particular, the leapfrog property) 
and so its analysis does not carry over to doubling spaces. In fact, the question of light low- 
stretch spanners for doubling metric spaces has resisted multiple attempts at resolution, and the 
problem of proving (or disproving) their existence has remained a central open problem in the 
study of spanners. The best lightness bound known for spanners in these metrics was H(logn) 
[Smi09l [ES13] . Our contribution is in proving the following theorem: 

Theorem 1.1. Given a graph G representing a metric S, there exists a (1 -|- e)-stretch spanner R 
of G with weight We ■ tc(MST(G)) where Wd = (ddim and 0 < e < ^. Spanner R can 
be constructed in time (ddim/e)‘^(‘^'^™)nlog^ n. 

Hence we resolve the question in the affirmative. We first prove that graphs which admit 
spanning trees that are everywhere sparse also admit light low-stretch spanners (Theorem 13.41 in 
Section [3|) . We then show that doubling spaces can be decomposed into sets that are everywhere 
sparse, and that the light spanners on these sparse sets can be joined together into a light low-stretch 
spanner for the original set (Section [4]). 


Related work. Light spanners are known for only a handful of settings, 
graphs [ADD+93[ lAGG+96 


unit disk graphs |KPX08] . 


These include planar 
graphs which are snowflakes of metrics 


|GS14j . and graphs of bounded pathwidth |GH12j and bounded genus [DHMlOj . 


Application to traveling salesman. A polynomial-time approximation scheme (PTAS) for the 
Euclidean traveling salesman problem (TSP) was presented by Arora [Aro98] (for all d-dimensional 
space) and Mitchell |Mit99| (for the Euclidean plane). These celebrated results were further im¬ 
proved upon by Rao and Smith [RS98| . who noted that the solution tour may be restricted to the 
edges of a light low-weight spanner for the sest. They used this observation to improve the runtime 
of the Euclidean PTAS. 
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Building on these frameworks, Bartal et al. [BGK12j showed that in metric spaces a (1 + e)- 
approximate tour can be computed in time roughly They also noted that whatever 

the minimum value of Wd may be, a tour can be computed in time log^ iog«- 77 . 

Immediate from the work of Bartal and Gottlieb |BG13] is a better runtime of 

plus the time required to compute the spanner. As a consequence of Theorem 11.11 we conclude 

that doubling spaces admit a PTAS with runtime 'n + (ddim/e)‘^(^'^‘™)relog^n. 

This is a dramatic improvement over what was previously known0 


1.1 Preliminaries and notation 

Graph properties. Let the weight of an edge e (rc(e)) be its length. Let the weight of an edge set 
E be the sum of the weights of its edges, w{E) = Similarly, the weight of any graph 

G = {Vg,Eg) is the weight of its edge-set, w{G) = w{Eg)- It follows that the weight of a path is 
its path length. For a partial graph R = {Vji,Eii) C G (meaning Vr C Vg, Er C Eg), diamG(R) 
is the diameter of R under the distance function cIg; it is the maximum distance in G between a 
pair of vertices also in R. Then diamR(i?) is just the length of the longest path in R. 

Let B{u,r) C Vg refer to the vertices of Vg contained in the closed ball centered at tt € Vg 
with radius r. B*{u,r) C Eg is the edge set of the complete graph on B{u,r). A graph G is said 
to be s-sparse if for every radius r and vertex v £ V, the weight of B*{v,r) n Eg is at most sr. 
Let A{u,ri,r 2 ) be the annulus that includes all points at distance from u in the range [ri,r 2 ], and 
A*{u,ri,r 2 ) the edge-set on the complete graph on A{u,ri,r 2 ). Points within distance n of u are 
within the hollow of the annulus, and not part of it. 

Doubling dimension. For a point set S, let A = A (S') be the smallest number such that every 
ball in S can be covered by A balls of half the radius, where all balls are centered at points of S. 
Then A is the doubling constant of S, and the doubling dimension of S is ddim = ddim(S) = log 2 A. 
The following lemma states the well-known packing property of doubling spaces (see for example 

[KL04]). 

Lemma 1.2. If S is a metric space and CCS has minimum inter-point distance b, then \C\ = 

Point hierarchies. Similar to what was described in [GGN061 IKL04j . a subset of points X CY 
is an r-net of Y if it satisfies the following properties: 

(i) Packing: For every x,y £ X, d{x, y) > r. 

(ii) Govering: Every point y £Y is strictly within distance r of some point x £ X: d{x,y) < r. 

The previous conditions require that the points of X be spaced out, yet nevertheless cover all points 
of T. A point in X covering a point in Y is called a parent of the covered point; this definition allows 
for a point to have multiple parents, but we shall assign to the point a single parent arbitrarily. 

^The result attributed to |BGK12 | is found in Theorem 4.1 of the first arxiv version. The focus of |BG13| was a 
fmear-time approximation scheme for Euclidean spaces under the word-RAM computational model, but aside from 
the time required to construct a hierarchy and spanner, as well as a MST for small areas, the entire construction 
holds for metric spaces. A theorem to this effect will be appended to a future version of |BG13| . and is beyond the 
scope of this paper. 
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Two net-points x, y of some r-net are called c-neighbors if d{x, y) < cr. By the packing property of 
doubling spaces, a net-point x can have c-neighbors. 

A hierarchy for a point set S is composed of nested levels of r-nets, where each level i of the 
hierarchy is a 2*-net of the level i — 1 beneath it. Note that the distance from an Tlevel point 
to all its descendants is less than ~ 2 • 2*. By scaling, we will assume throughout this 

paper that the minimum inter-point distance in all sets is at least 1. Then we have that top 
level H := [log 2 diam(iS')] contains a single point, and the 0-level contains all of S. (Below, we 
will find it convenient to allow for levels i < 0, and each of these levels contains all of S as well. 
The bottom level is L) A full hierarchy for S constructed in time 2‘^('^'^‘“^relog A, where A is the 
aspect ratio of S, the ratio between the largest and smallest inter-point distance in S [KLOd] . A 
hierarchy can also be constructed in time 2‘^*^‘^'^‘“)relogn, and in this case some isolated net-points 
are represented implicitly [HM061 [CG06| . We can also augment the hierarchy with c-neighbor lists 
for all hierarchical levels and points: It is known that there are only ^^(ddim)^ neighbor pairs, and 
that given the hierarchy they may all be discovered in (P{Ad\m)^ time and space. 

An r-semi-net is a net that satisfies the covering property but not the packing property. A 
semi-hierarchy is composed of levels of nested semi-nets. 

Net-Respecting spanning trees. We will say that a spanning tree T is net respecting (NR) 
relative to a given hierarchy if for every edge e in T, both of its endpoints are f-level net-points 
for i satisfying 12 • 2* < w{e) < 24 • 2*"*“^. (Recall from above that we may assume the existence of 
levels below i = 0.) The following lemma is adapted from Lemma 1.6 in [GKK13] . 

Lemma 1.3. Every spanning tree T can be converted into a net-respecting spanning tree T' that 
visits all points visited by T, such that w{T') < 5w{T). 

Proof. For every edge e = {x, y) in T do the following. Let x', y' be Tlevel net-point ancestors of x, y 
respectively, for the lowest i satisfying 12 • 2* < d{x',y') < 24 • 2*+^. Replace edge e with long edge 
e' = (x',y'), and also add the short edges {x,x') and {y,y'). The short edges have weight at most 
2 • 2*, while by the triangle inequality w{e) > w{e') — 4-2* >8-2*. While edge e' is net-respecting, 
edges (x, x') and (y, y') may not be. The short edges are themselves replaced by the procedure 
above. This leads to a series of edge replacements, and if the replacement adds an edge already in 
the set, we will remove the duplicate. The total weight of all short edges added by the recursive 
procedure is a geometric seriers summing to less than 8 • 2* < w{e'). As the ratio of the long edge 
to the original edge is at most |, the total weight of all edges is 2w{e'){l + |) = 5w(e'). □ 

Let MST^^ denote the net-respecting minimum spanning tree. The following lemma related a 
local minimum spanning tree to a global net-respecting minimum spanning true. It is immediate 
from Lemma 1.11 in |GKK13| . which made a similar claim concerning TSP tours. 

Lemma 1.4. Let S be a point set equipped with a hierarchy, and S' C S. Let u he any point of S, 
and r > 0 any value. Then 

(i) w{MST^^{S)nB*{u,r)) < 14(t(;(MST(R(u, r))) 

(ii) w{MST{B{u,r))) < 2w{MST^^{S) D B*{u,4r)) + 

We recall that the removal of points from a set can increase the weight of its minimum spanning 
tree by at most a factor of 2. Let S' consist of B{u, 2*) along with all hierarchical ancestors of these 
points up to level i. Then S' C B{u, 3 ■ 2*). We conclude that 

l^w{MST^^{S') n B*{u, r)) < w{MST{S')) < 4w{MST^^{S) n B*{u, 12r)) + (1) 
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Finally, we note that for any S, t(;(MST(S')) < IS"! diam(S'). Tighter bounds are in fact known 
for doubling spaces [Aro98t [Tain4[ [Smiin| . but are not necessary for our results. 


2 Background: Hierarchical spanners 


Before presenting our results, we review the basic hierarchical spanner introduced by Gao et ah 
[GGN06] . along with some small modifications suited for our purposes. This spanner provides a 
theoretical framework for our approach, although our ultimate construction will be significantly 
more involved. The results presented in this section are mostly well-known. 

Throughout this paper, we will make the simplifying assumption that G is a complete graph. 
A complete hierarchical spanner R for a graph G on 5 is derived from the hierarchy of S. For some 
constant c, for each levels i we add an edge between every pair of Tlevel net-points with inter-point 
distance at most c2*. More precisely, we add to Er an edge between the two vertices in Vr that 
represent these points. Recall that the top level of the hierarchy is H, and the minimum inter-point 
distance in S' in 1. For this construction we will assume that the bottom level of the hierarchy is 
L : [log(l/c)J. Let M = H — L-\-l he the number of levels in the hierarchy for S. To construct the 
complete hierarchical spanner, it suffices to construct a hierarchy for S and compute all relevant 
c-neighbors, in total time 2‘^(^‘^™)nlogn -|- [GROSbl IGROSa] . 


Lemma 2.1. Any spanner R for graph G on S which satisfies that all i-level net-point c-neighbor 
pairs have 0-stretch (for all i), is a {1-\- ^)-stretch spanner for all of G, when c > 24. 


Proof. Consider any two net-points x,y in the bottom level L of the hierarchy of S. Let x',y' be 
their respective z-level ancestors, where i is the minimum value for which x',y' are c-neighbors 
We have dc{x',y') < c2L To derive a lower-bound on dcix',y'), let x",y" be the {i — 1)- 
level ancestors, and by assumption dcix",y'') > c2*“^. Applying the triangle inequality, we 
have dcix',y') > —dcix'^x") -\- dcix”.,y") — dciy'',y’) > c2*“^ — 2-2* = (c — 4)2*“^. Since 


x',y' are c-neighbors, dR{x',y') = dcix'^y'). Applying the triangle inequality we have < 

dji{x,x')+dR{x',y')+dR{y',y) dgix',y')+2-2^+^ (c-4)2"~ _ e-|-4 _ i 16 < i i ^ n 

—da{x,x')+da(x',y')—dG{y',y) dQ{x',y')—2-2'‘+^ (c—4)2®“i —2®+^ c—12 ' c—12 — ' c ' 


The stretch of this spanner is arbitrarily low, but it has poor lightness bounds: In the trivial 
case where the n points of S reside on the line at intervals of distance 1, it is easy to see that 
w{Er) = 0(nlogn) = 0(logn• u)(MST(iS')). We may consider being more stingy with added edges 
by using the greedy algorithm of [DN97] (see also [GLN02| i. and we call this construction the greedy 
hierarchical spanner: We build the spanner as before, but first sort all cOfddim)^ edges in increasing 
order. We consider each edge in turn, and before adding an edge between c-neighbors, we check 
to see if the stretch between them on the current partial spanner exceeds some threshold - for 
example, it is at most 1-1-6 for some constant 0 < 6 < 1 - and only add the edge if the condition is 
met. In Lemma 12.21 below, we show that the greedy hierarchical spanner has 1 -|- 0{^ -p b) stretch. 


Lemma 2.2. A spanner R for graph G on S, for which all i-level net-point c-neighbor pairs have 
(1 -|- b)-stretch (for all i) for any constant 6 > 0, is a (1 ^ 6b)-stretch spanner for all of G, 

when c > 24. 


Proof. The proof is similar to the proof of Lemma 12.11 Consider any two net-points x, y in the 
bottom level L of the hierarchy of S. Let x', y' be their respective i-level ancestors, where i is the 
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minimum value for which x',y' are c-neighbors: dcix',y') < c 2 ®. As above, dcix',y') > (c —4)2* 
Now noting that by construction dji{x', y') < (1 + b)dG{x', y'), and applying the triangle inequality, 

< dji{x,x')+dji{x',y')+dR{y',y) {l+b)dGix',y')+2-2'^+^ (l+b)(e-4)2*-l+2*+^ _ e+4 i 

dG{x,y) — —dax,x'+dG{x',y')—dG{y',y) dcix',y')—2-2’‘+^ (c—4)2* 2®+^ c—12 ' 

+ f + □ 

In closing this section, we note that the stretch bounds of Lemmata 12.11 and 12.21 also hold for 
net-points of semi-hierarchies. Indeed, the proofs use only the covering property of nets, and not 
their packing property. 

3 A light spanner for graphs with sparse spanning trees 

In this section, we show that graphs with sparse spanning trees admit light low-stretch spanners. 
Note that we do not require the graph itself to be sparse, merely that it admit a sparse spanning 
tree. In Section EH we will show that we can build a light spanner for the union of two low-stretch 
paths. In Section [3.21 we will build on this result to build a light spanner for all graphs with sparse 
spanning trees. 

3.1 A light spanner for pairs of close low-stretch paths 

In this section we show that given a close pair P,QcGof low-stretch paths, we can compute a 
light low-stretch spanner for their union. We first need to define nets, hierarchies and bipartite 
spanners for paths. 

Path nets and hierarchies. Let P C G be a path consisting of some vertices and edges in G. 
An r-path-net for P is an r-net for P under the path distance function dp. A path-hierarchy for P 
is a full hierarchy under dp, with net points in levels i = [log 2 w{P)~\ and lower. 

It is easy to see that an r-path-net of P is an r-semi-net for the point of P under do- We 
observe that the distance function dp is non-contractive with respect to do] dG{u,v) < dp{u,v) for 
all u,v ^ V. Hence, if a point is covered in the path-net, it is covered under dG as well. However, 
the packing property may not hold, so the path-net is only a semi-net for P under dG- It follows 
as well that the path-hierarchy for P is a semi-hierarchy for P under dG- 

Bipartite path spanners. Let G = {Vg,Eg) be a complete graph, and let P = {Vp,Ep),Q = 
(Vg, Eq) be two paths in G with stretch at most (1 -|- hi) for some 1 < 6 i < c. Let each path be 
equipped with a net-hierarchy. A complete bipartite hierarchical spanner R for P U Q contains all 
vertices of Vp U Vq and path edges Ep U Eq as well as the set of all edges E C Eg connecting 
vertices of Vp and Vq . A greedy biparite hierarchical spanner R for P Q contains all vertices of 
Vp U Vq and path edges Pp U Eq, as well as an edge-set E' C E C Eq- We sort the edges of E in 
increasing order, and consider each edge in turn. An edge of E is added to E' only if the stretch 
between its endpoint on the current partial spanner is large, say more than (1 -|- 5i -|- 62 ) for some 
0 < 62 < 1 — bi. We show that the greedy bipartite spanner has favorable properties: 

Lemma 3.1. Let P,Q C G be a pair of (1 -|- bi)-stretch paths whose distance from each other is 
not greater than c ■ min{r(;(P), (for 0 < 61 < 1 and any c > 24). Let A be the aspect ratio of 

P L) Q. A greedy bipartite hierarchical spanner R with parameters c and 0 < 62 < 1 ~ for paths 
P,Q can be constructed in time 0{cs{\Vp\ + |Vq|) log(cA)), and satisfies the following properties: 

(i) Stretch: dp{p, ( 7 ) < (1 ^ -k 6(61 -|- b 2 ))dG{p, q) for all p & P and q ^ Q. 
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(a) Weight: w{E') = where E' = Er — {Ep U Eq) is the set of new 

spanner edges not in P or Q. 

Proof, (i) This follows immediate from Lemma 12.21 

(ii) The proof is by a charging argument. We assume without loss of generality that w{P) < w{Q). 
Take the maximum level i for which S has an i-level edge, and let the endpoints of the edge be 
p ^ P and q G Q. Consider any point p G P satisfying dp{p,p) < r for r = ^dG{p,q), and we will 
show that p cannot have any edge incident to Q in Ep: 

First note that the proximity of p to p implies that p may be a j-level net-point only for values 
j < logr. Now suppose by way of contraction that p G P has an edge to some q G Q, and then by 
construction in must be that dciP, q) < c2^ < cr. Let R' be the partial spanner before the addition 
of the i-level edges. Then the stretch guarantees of the paths and spanner, along with applica¬ 
tions of the triangle inequality, imply that ^ 

r-l-er-l-(l-l-bi)[dg {p,q)+da{p,p)+dG (p,g)] ^ {2+bi){r+cr)+ll+bi)dG{p,q) _ i , . , {2+bi){r+cr) , , , Acr _ 

<iG(p,q) — dG(p,q) “T IT" dGip,q) ^ dglp,!}) 

1 + bi + b 2 . This implies that in R' the stretch from p to g is not sufficient to add to E' an edge 
connecting them - a contradiction. 

It follows that for any i-level edge e G E' connecting p G P to q G Q, the entire path of P within 
distance ^w{e) of p under dp has no lower level edges incident upon it. We charge the Llevel edge 
to this segment of the path P. 

The statement now follows by noting that p may have at most 3cs i-level edges incident upon 
it (all of which will be charged to the same path segment of P): This is because all relevant i-level 
path-net points of Q within distance c2* of p. Since the stretch in Q is less than 2, a path connecting 
all these points mush lie fully within distance 3c2* of p, and by the the sparsity guarantee its length 
is less than 3cs2*. There can be at most 3cs i-level path-net points on a path of this length. 
Runtime. At each of 0(log(cA)) levels i, we must locate for each i-level path-net point p G P 
its i-level path-net c-neighbors in Q. As shown above, there are most 3cs such neighbors. Further, 
if two points p G P and q G Q are Tlevel c-neighbors, then so are their respective parents p' G P 
and q' G Q: dcip'^q') < dcip^p') + dcip^q) + dciq.q') < 2 • 2*+^ c2* = (f + 2)2* -|- 1 < c2*+b 

So to compute all c-neighbors, it suffices to iterate down the hierarchy, maintaining for each i-level 
net-point p G P a, list of all c-neighbors in Q. The lists for the children of p can be found by 
considering all children of p’s c-neighbors, for a total of 0(cs) candidates. □ 

3.2 Extension to graphs with sparse spanning trees 

The results in the previous section apply only to pairs of low-stretch paths. We will use them to 
obtain similar results for all graphs with sparse spanning trees. Our plan is as follows: We first 
show how to decompose a spanning tree into paths without a stretch guarantee iLemma 13.21) . We 
then show how to replace any path with a light set of low-stretch paths (Lemma 13.3p . Finally, we 
build the bipartite spanner of the previous section on all pairs of low-stretch paths, and show that 
this gives a light low-stretch spanner tree for the full set (Theorem 13.4p . To do this, we will need 
to demonstrate that the union of all the low-stretch paths is still relatively sparse, and also that 
the union of the path-nets for these paths can serve as a semi-hierarchy for the entire space. 

Lemma 3.2. Given an s-sparse spanning tree T = (Vr, Ep) C G of n nodes with aspect ratio A, 
T may be decomposed in time 0(s|VT|logA) into a set Q of paths, with the following properties: 

(i) Diameter: At least one path Pi G Q has length diam 7 ’(r). 
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(ii) Proximity to path-net points: Every vertex in Vq is within distance b (under dr) of some path 
Pj ^ Q of length w{Pj) > b, for all values 0 < b < diam'r(T') . 

Proof. The decomposition procedure removes from the spanning tree the longest path - of length 
diam(P) - and places the path in the collection Q. That is, the edges of the path are removed 
from Et, and then all vertices with no edges incident upon them are removed from Vt- After the 
removal of the longest path, a number of disjoint subtrees may remain, and each is decomposed 
recursively until Vt is empty. This completes the description of the decomposition procedure. 

We note trivially that the distance under dx from a removed path to all vertices in one of the 
remaining spanning trees is not greater than the graph diameter of that spanning tree. 

(i) The first item follows by construction. 

(ii) Take any vertex and the paths that were removed from the spanning tree containing the vertex 
at each iteration. The paths are necessarily in decreasing order of length. Consider the first path 
of length at most b. Then the diameter of the subtree that contained this path before its removal 
was at most 6, hence the the vertex is within distance b of the next path, which has length greater 
than b. 

Runtime. First note trivially that longest path of a tree can be found in 0{n) time. Now consider 
any subtree formed by the removal of a path P of length w{P). After 0(s) further decomposition 
steps, the longest path in any remaining subtrees cannot be longer than . This is because each 
decomposition step removes the longest path of the current subtree, and leaves behind it a group of 
subtrees all within distance w(P) of P. Then the fact that the original subtree is s-sparse implies 
that it may contain at most 0(s) paths of length greater than It follows that the total runtime 
of the decomposition procedure is 0{s\Vp\ log A). □ 

It remains to replace each of the paths in Q with a set of low-stretch paths. Consider the follow¬ 
ing procedure (motivated by the greedy algorithm) which takes a path P and value c, and replaces 
P with a set of low-stretch paths V. First consider the endpoint p,q £ P. If the endpoints are 
close together - dp{p,q) < idiamG'(P) - find a point r satisfying dG{p,r),dc{q,r) > |diamG(P), 
segment P at r, and continue the procedure below separately on the two smaller paths. Otherwise, 
continue the procedure on P itself; 

Iteratively create a path-net for P, at each iteration i = L, L 1,... promoting some of the 
(i — l)-level path-net points to also be i-level path-net points; We begin by promoting the first 
point in P - call this p - and proceeding down the path promote in turn every (i — l)-level point 
at path distance 2* or greater from the previously promoted point. 

Then beginning again at p, we search down the path for the associate of p, the first i-level 
path-net point q satisfying dG{p,q) > c2*. Now, if the stretch is greater than 1 -|- , 

remove the partial path connecting p, q, and replace it by a directed edge. The partial path is then 
segmented into smaller paths at its (i —2)-level path-net points, and these small paths are all added 
to V. If the partial path is removed, the procedure continues with point q, searching down the 
path for its associate. If the partial path is not removed, then the procedure takes the first i-level 
path-net point r following p, and searches down the path for r’s associate. Iteration i terminates 
upon reaching an i-level path-net point r with no associate. In this case, we locate the final i-level 
path-net point in P ~ call this q - and remove the partial path connecting p and q, replacing it 
by a direct edge. The removed partial path, as well as the path beyond q, are decomposed into 
smaller paths at their (i — 2)-level path-net points, and added to P. The procedure terminates 
upon reaching a level i for which the endpoint p has no associate. 
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Lemma [3.31 below shows that the final path set V has favorable properties, which we will use in 
the construction of the spanner. 

Lemma 3.3. Let P = {Vp,Ep) be an s-sparse path in G with arbitrary stretch and aspect ratio A. 
For any constant c > 24, the collection V may be computed in time 0(cs|4^| log A), and possesses 
the following properties: 

(i) Vertex cover: The union of vertices in all sets Pi ^ V is exactly Vp (ViVp^ = Vp), and 
EJVpJ<3|Vp|. 

(a) Stretch: Each path in V has stretch at most (1 + ^). 

(Hi) Diameter: At least one path Pi has length ^ or greater. 

(iv) Proximity to path-net points: Every vertex v G Vp of is within do distance 16(c+ 1)6 of some 

path Pj & P of weight w{Pj) > b, where b is any value satisfying 0 < 6 < _ 

(v) Path sparsity: Each path in V is 2>s-sparse. 

(vi) Graph sparsify and weight: The union of all paths Pi gP form a graph that is both 3s(3cs+l)- 
sparse and has total weight at most 3s(3c(c + l)s + 1) • w{P). 

Proof, (i) This covering immediate from the construction, and the number of vertices follows by 
noting that a removed path may duplicate the path endpoints, but removes at least a single vertex 
from the original path. 

(ii) Let p, q be any pair of z-level path-net points in a path P' G P, and we show that the path 
stretch of p, q is less than 1 + ^• Then the item will follow from Lemma 12.21 

Since p, q are found on the same path, it must be that all (i -|- l)-level associate pairs r, t 
straddling the path segment from p to q have path stretch less than 1 -|- , as otherwise the 

path connecting r and t would have been removed and segmented along its (i — l)-level path-net 
points, and p, q could not be in the same path segment in P. Recall by construction the dc- 
distance from r to t is in the range [c2*"*“^, (c-|- 1)2*'’'^), and since all paths of P are 3s-sparse 
(item (iv) of the lemma), the dpi distance from r to t is at most 3s • (c -|- 1)2*"*'^. Were the 
pair p, q to have path stretch at least 1 + ^) then r, t would have path stretch at least ~ 

dpi {r,p}+dpi {r,t)+dp, (t,p) . dc(rp)+{l+\)dG(r,t)+dpi (t,p) dG{r,t)+^)dG(r,t) ^ 3s-(c-|-l)2"+W2Vc _ .. 

do {'I'd) — dcVd) — dG{r,t) — 3s-(c+l)2®+^ ' 

3c(c+i)s ’ which is a contradiction. 

(iii) The endpoints p,q G P were chosen so that diamp(P) > gy construction, the final 

path remaining after all iterations retains p as an endpoint. On the other end, at each iteration 
i a small end-segment of the path may be segmented, but these account for less than a fraction 
c E£o ~ total path length. So the diameter is at least ^ > diamG(P) ^ 

Finally, the s-sparsity of path P directly implies that diamp(P) < sdiamc'(P). 

(iv) For every added edge connecting i-level path-net points p, q, the removed path was all within 
distance c2* of p (under da) of the new edge connecting p, q, which have length in the range 
dcip, q) G [c2*, (c-|-l)2*). Now the removed path may itself have had paths removed from it at earlier 
stages of the iteration, but all these paths must be within distance E}=o < c2^^^ < 2dG{p, q) of 
the new edge. Also note that when this path was removed, it was decomposed into smaller paths of 















path-length at least 2*“^ = Hence, the length of the edges feature jumps of at most 8 (c -|- 1), 
and so every vertex is within distance 16(c -|- 1)6 of an edge of length at least 6 . 

(v) Take an edge added by the procedure between points p, q. Then in the original path P there 

must have been a heavier path connecting p, q, and by the arguments in item (iv) this path was 
fully within distance dcip, q) of both p and q. Any ball containing the new edge must be of radius 
at least and so it contains the entire (removed) old path within 3 times its radius. So the 

weight of the smaller ball is less than the weight of the larger larger ball, and the sparsity is 3s. 

(vi) For the claim of weight: We have already noted in item (iii) that the edges added when 

removing an end-segment sum to much less than the weight of the path. For the other edges, the 
proof proceeds by a charging argument. For an edge e added between points p, q in iteration i, 
charge its weight to the edges of the removed path. Each edge e' in the removed partial path Q 
receives change and since Q must have weight at least u;(e)(l -|- 3 c(e+i)s )i ^he charge to e' 

is at most ^ —. Now, if e' itself was added in an earlier iteration j, then the new charge to e' 

3c(c+l)s 

due to the addition of e is itself passed down to earlier removed paths, until all charges lie solely on 
edges present in the original path P. Let e be an edge of P, and is follows that the sum of charges 
placed on e (included the cost of e itself) is less than w{e) - 3 c(c+i)s ) * < (3c(c+l)s+l)u;(e). 

For the claim of sparsity: It follows from the proof of items (iii) and (iv) that all edges of V 
inside an r-radius ball charge to edges of the original path P within distance 3r of the ball’s center. 
As P is s-sparse, the total weight in the small ball is 3s(3c(c -|- l)s -|- 1), and the total weight of all 
edges in V is 3s(3c(c -|- l)s -|- 1) • w{P). 

Runtime. There are 0(log A) iterations, and at each one we can make a single pass on the path, 
promoting net-points and recording their order in a list. Then we must discover for each i-level 
path-net point p its c-neighbors. But these are all within path distance 0(cs2*) of p, so it suffices 
to inspect only the next 0(cs) net-points in the list. The runtime follows. □ 

Joining together Lemmata l3.1l[T^ and l3.3l we can prove that spaces with sparse spanning trees 
admit light low-stretch spanners: 

Theorem 3.4. Let T = {Vt,Pt) be an s-sparse spanning tree of a complete graph G with aspect 
ratio A. Then in time (s/e)*^^^^ • \Vt\ log A we can construct for G a spanner R with 

(i) Weight: VF, • w{T) for W, = (s/e)^W. 

(a) Stretch: (l-|-e). 

Proof. The construction is straightforward: Decompose the spanning tree T into paths using the 
decomposition procedure of Lemma 13.21 and then replace each path with a set of low-stretch paths 
using the replacement procedure of Lemma 13.31 with parameter ci = 3 • 2^/e. For every pair of 
paths in the new set, build the greedy bipartite spanner of Lemma 13.11 with parameters C 2 = 
max{16ci -|- 18,8s} • 2®/e and 62 = e/12. The new graph is G'. Recall that the bipartite spanner 
will only add edges if the path pair is sufficiently close. 

(i) We first must prove sparsity: Tree T is s-sparse, and the sparsity condition is not affected by the 
decomposition. Take any path P the decomposition, and P is replaced by a set of paths connecting 
the vertices of P with sparsity O(cfs) (Lemma 13.3l vi)). Hence, before the addition of the bipartite 
edges G' possessed sparsity s' = 0{c\s‘^). 

Returning to proof of weight, it follows from Lemma l3.3l vi) that the weight of the paths of G' 
before the construction of the bipartite spanners was 3s(3ci(ci -|-l)s-|-l) •rc(T). A bipartite spanner 
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with parameter C 2 is built for each pair of paths, and we will charge the cost of the new edges to 
the shorter path. By sparsity, any path P is within distance C 2 • w{P) of 0{c2s') = 0 (cfc 2 S^) 
longer paths. Charge P for the additional edges added by the bipartite spanner to each of these 
paths. By Lemma [3B] this yields a charge of 0 (c|s/ 62 ) • w{P) per spanner, for a total charge of 
0(cf€ 25 ^/ 62 ) ■ w{P). It follows that w{G') = 0(cf 025 ^/ 62 ) • w{T). 

(ii) By Lemma 18.31 the intra-path stretch in all paths produced by the replacement procedure is 
1 By construction of the bipartite spanner, the stretch among all i-level path-net points that 

are C 2 -neighbors is max{l -|- 62 ) 1 + — }• 

We require a semi-hierarchy: Below we show that each vertex v in G is within distance (16ci -|- 
18)2* of some Llevel path-net point for all levels i < [diam'r(r)/ 8 s]. It follows that if we promote 
all i level path-net points (for each i) to be j-level points for j = i + logmax{16ci -|- 18,8s}, 
then each vertex v is within distance at most 2^ of some j-level point, and this hold for all levels 
j < [diam 7 ’(T)] = H. Hence we have a valid semi-hierarchy for all points. Substituting j into the 
stretch bound above, we have that the the stretch among all j-level semi-hierarchy net-points that 
are max{i 6 ci+i 8 sg} ~ii 6 ighbors is max{l-|- 62 , l + Substituting in the values of ci, 62 chosen above, 
we have that the the stretch among all j-level semi-hierarchy net-points that are (64/e)-neighbors 
is 1 + ^. Then by Lemma [221 the resulting spanner has stretch 1 -|- -|- 6 ^ = 1 + e. 

It remains only to verify the above claim, that each vertex in G' is within distance (16ci -|- 18)2* 
of some (not yet promoted) i-level path-net point: After the decomposition procedure, for any 
value 0 < 6 < diamr(T) each vertex v € Vq is within distance b of some vertex w a Pi for which 
w{Pi) P b (Lemma I3.2l ii)i. After the replacement procedure, w is within distance 16(ci -|- l)b' of 
some path of length at least b' where Q <b' < ^ p = have that v is within 

distance (16ci -|- 17)6 of some path of length at least 6 , where 0 < 6 < ^ path of this 

length possesses a path-net point at level i = [log(6/2)] or greater. It follows that v is within 
distance (16ci -|- 18)2* of some Llevel path-net point for all levels i < [diam'r(T)/ 8 s]. 

Runtime. By Lemmata 13.21 and 13.31 the decomposition and replacement steps can be done in 
time 0{c2s\Vt\ log A). To compute each bipartite spanners, for each path P in G' we must have all 
paths of length w{P) or greater within distance C 2 ■ w{P) of P. Recall from the proof of (i) above 
that the sparsity of G' before the addition of bipartite spanner edges is s' = {s/e)^^^\ 

After the initial decomposition step, assign a path-hierarchy to each path, and compute all 
i-level 2 c 2 -neighbors among all paths and pairs of paths, for all i. This can all be done in time 
(s/e)*^^^^?! log A. (The construction and analysis is identical to that of Lemma 13.2B For each 
path e created by the replacement step, we compute its distance to all Tlevel path-net points in 
the original edge within distance C 22 *. Since the current graph is s'-sparse and has aspect ratio 
A this can be done in time (s/e)*^^^)/! log A over all paths, if we precomputed for each vertex its 
closest Tlevel path-net point for all i. For each i-level path-net point, we maintain a list of the 
replacement path which discovered this point, sorted in order of path length. Note that we may 
sort all replacement paths once, and then compute close path-net point for each one in turn. 

For each path P and each i < log(rc(e)), we inspect the i-level path-net points within distance 
C 22 *, and all z-level C 2 -neighbors of those path net points, and their listed paths of length greater 
than w{P), and this suffices to discover for P all longer edges within distance C 2 -w{P). The runtime 
follows from the sparsity and aspect ratio of the current graph. □ 
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4 A light spanner for general metric graphs 

In Section [3] we showed how to construct a light spanner for spaces that have sparse spanning trees. 
In this section we complete the proof of Theorem II.II bv showing how to decompose general metric 
graphs of low doubling dimension into graphs with sparse spanning trees. Then the light spanners 
for the decomposed sparse spaces can all be joined into a single light spanner for the original metric 
space. Our approach uses a technique developed for computing near-optimal traveling salesman 
tours in polynomial time [BGK121 IBG13| . although our setting is much less restrictive, and so the 
problem of finding a good decomposition is simpler. 

For some graph G equipped with a hierarchy, define F{u,i) for i-level vertex u to include 
all points of B{u,i) n G, as well as all their hierarchical ancestors up to level i. The following 
preliminary lemma shows that we can spin off a (slightly) dense area of the graph: 

Lemma 4.1. Let G be a graph with aspect ratio A, equipped with a hierarchy. For some fixed f = 
(ddim ^ ^g ^^g iQyjQgi level for which there exists some u satisfying MST^^{F{u, i)) > 

/2*. Then G may be segmented into intersecting subgraphs G', D G G with the following properties, 
for any value c> 2: 

(i) Weight in G’: w{MST^^{F{u,i)nG')) < ^2\ while wiMST^^{G))-w{MST^^iG')) > 

(a) Sparsity of D: MST(ZI) is f-sparse, while D has diameter 0{c) ■ 2*. 

(Hi) Neighbor proximity: For every c-neighbor point pair in the hierarchy of G, the pair is found 
together in A or B or in both. 

Proof. Subset D includes all points in B{u, (13 -|- c) • 2*), along with all their net-points up to level 
i. Subset G' includes all points of G — B{u, 13-2*), as well as all z-level (and higher level) net-points 
in B{u, 13 • 2*). Below, we will also add some more points to A. 

(ii) Recall that by assumption each {i — l)-level point v satisfies MST^^{F{v,i — 1)) < /2*“^. A 

spanning tree for any ball of radius 2* < r < (13 -t- c) • 2* in D can be formed by covering that ball 
with the spanning trees of sets F{v,i — 1), and then connecting the centers of these sets 

at an additional cost of for a total weight of By the optimally of MST{D), 

its weight inside the given ball cannot be greater than this. A similar result holds for values of r 
less than 2*, covering with sets F{v,k) for k < i — 1. 

(iii) Glearly c-neighbor pairs of levels i or higher are both found in G'. For levels k < i, ii one of 
the points is not found in G', then by construction it is found in B{u, 13 • 2*), and so the second 
point is found in B{u, 13 • 2* -|- c2^) C B{u, (13 -|- c) • 2*) = D. 

(i) We first discuss the weight of MST'^^(G) in the area around u: A consequence of Equation ([1]) 
is that B*{u, 12 • 2*) had intersected edges of MST^'^(G) weighing more than ^2*. At the same 
time, the larger ball i?(ri, 13 • 2®) admitted a minimum spanning tree of weight (by 

arguments identical to the proof of item (ii) above). So by Lemma ll.dl if B{u, 13 • 2*) intersected 
edges of MST^^{G) of weight 2‘=’('^^'”)/2h 

We now add some more points to G': Set j = i — a log ddim for some absolute constant 
a > 2 to be specified below. Add to G' all j-level points in B{u,13 ■ 2*). Further, let radius 
r G [12 • 2® + 72 • 2fi 13 • 2® - 72 • 2^] be such that u;(MST^'^(G) n A*{u, r - 72 ■ 2fir + 72 • 2^) < 
^w{MST^^{G))r\B* {u, r—72-2^)) - that is, the weight of MST^^(G) inside the annulus is a fraction 
of that inside the hollow. (Such a value r must in fact exist: There are 0(2®/2'^) = 0(2“ ddim) 
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non-intersecting annuli. If each contained edges of MST'^'^(G) weighing a fraction at least | of the 
edge-weight within their hollow, the total edge-weight within B{u, 13-2*) would be 
for a sufficiently large choice of a, this exceed the upper-bound proved above.) Add to G' all points 
of the smaller annulus A(u,r — 24 • 2-^, r -|- 24 • 2-^), as well as all points outside outer edge of the 
annulus - in short, all points outside B{u,r — 24 ■ 2^). 

To prove the first part of item (i), note that F{u, i)r\G' contains only points of level j or higher, 
and these can connected by an MST of weight = ddim*^^^”"2* ^2*. The 

final inequality follows when / is chosen to be sufficiently large with respect to a. Then by Lemma 
o a connecting net-respecting MST has weight less than ^2*. 

To prove the second part of item (i), we will construct a net-respecting spanning graph for G' 
with smaller weight than MST^^(G). Recall that G and G' differ only on points within B{u,r — 
24-2®). We add to the spanning graph of G' all edges of MST'^'^(G) fully outside of B{u,r+ 24-2^) 
(and this is outside the smaller annulus A{u, r — 24-2^ ,r + 24-2^)). We also add to G' those edges 
with only a single endpoint inside R(n, r -|- 24 • 2®), and these net-respecting edges must be incident 
on the annulus A{u, r — 24 ■ 2^ ,r + 24 ■ 2^) if they are shorter than 48 • 2-^, or on j-level points in 
{u, r — 24-2^) if they are longer. We then claim that the points of G' inside the ball B{u,r + 24-2^) 
can be covered by a graph whose weight is lighter than that of MST(G) H B*{u, r + 24 ■ 2®) by at 
least ^2®, from which the second part of the item follows. 

First note that a graph connecting all j-level points in B{u, 13 • 2®) has weight much less than 
^2® for appropriately chosen / (as above in the proof of the first part of (i)), and so by Lemma 
Ol a net-respecting MST connecting only these points has weight much less than ^2®. Further, 
all points of the annulus A(u, r — 24 ■ 2^, r + 24 ■ 2^) can be connected to some j-level point at total 
cost at most w{MST^^{G) fl A*{u, r — 72 ■ 2^ ,r + 72 ■ 2^)) (which was chosen above to be at most 
ir(;(MST'^^(G)nR*(u,r —72-2-^))): The annulus cuts w{MST^^{G)) into disjoint paths. For each 
such path, if it is incident on a j-level (or higher level) point then we are done. Otherwise, consider 
the tail of this path, the final segment exiting the smaller annulus A(u, r — 24 • 2-^, r -|- 24 • 2-^). Since 
MST'^'^(G) is net-respecting, and the path touches only net-points of level less than j, the tail must 
be of length at least 24 • 2-^, since it cannot jump immediately to a point outside the larger annulus. 
Cut off this tail, and instead add to the path a much shorter tail of length at most 2 • 2^ connecting 
to the closest j-level point. So the total cost of the spanning graph for G' 0 B{u,r + 24 ■ 2^) is at 
most j^2® -|- jw{MST^^{G)r\B*{u, r — 72-2^)), while MST'^^(G) nR*(u, r-|-24 • 2-^ ) is of course at 
least u;(MST^^(G) n B*{u, r - 72 ■ 2^)), a difference of lw{MST^^{G) n B*{u, r - 72 • 2^)) - ^2® > 
lw{MST^^{G) n B*{u, 12 • 2^) - ^2® > f • {2® - ^2® = |2L □ 

By using the technique of Lemma [4.1l to repeatedly spin off relatively sparse areas of the graph, 
we can decompose the entire graph into relatively sparse areas, as in the following theorem: 

Theorem 4.2. Any graph G with aspect ratio A and equipped with a hierarchy, may be segmented 
into a set of subgraphs V (UiDi = G) with the following properties, for any c > 2 and for f as 
above: 

(i) Sparsity: For all Di G V, MST(A) is /-sparse. 

(ii) Weight: J2iw{MST{Di)) = c^^'^‘^^^w{MST^^{G)) 

(Hi) Neighbor Proximity: For every c-neighbor point pair in the hierarchy of G, the pair is found 
together in at least one set Di € V. 
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(iv) Point occurrence: Each point of G appears in at most c'^('idim) different sets in P. 

The segmentation can he computed in time c‘^(^'^‘“^relogn(logn + log A). 

Proof. The technique of Lemma l4. II requires knownledge of MST^^(G) (to compute the annulus), 
so we must first build this graph. We construct the full hierarchical spanner with parameter c. 

This can all be done in time log n, and the resulting graph has only edges 

[GROSbl IGROSaj . We then extract from the edges only those that are net-respecting, and build 
MST'^^(G) using this set, in total time log n. We will need to maintain MST^^(G) under 

a long series of point deletions, and this can be done in total time c‘^(‘^'^‘™)nlog^ n (see Theorem 6 
in [HdLTnij l. 

Beginning at the lowest hierarchical level L and iterating upwards, for each i-level net-point u 
we compute F{u, i) and construct MST^^(F(n, i)), applying the segmentation technique of Lemma 
lO if necessary. The subset D is added to V, and we iterate the procedure on subset G' until there 
are no more heavy neighborhoods, at which point the final subset G' is added to P as well. 

The sparsity of sets in P and neighbor proximity follow immediately from Lemma l4.1f iil and 
(hi), respectively. The bound on point occurences follows from the packing property and height of 
the hierarchy, and a corollary of this bound is that c‘^(^'^™^nlognlog A time suffices to compute 
all graphs MST^^{F{u,i)). For the weight, let Gi be the graph in the z-th iteration, and Wi = 
u;(MST^^(Gi+i)) - w{MST^^{Gi)). Glearly ^ • u;* < MST^^(G). By Lemma EKi) and (ii), 

^. w{MST{Di)) = = c®(^dim)y^^MST^^(G)). □ 

We can now complete the proof of Theorem I l.R 

Proof. Fix c = 64/e. Given graph G, we build a full hierarchy for the entire space, and add to 
the spanner R the edges of the complete hierarchical (1 -|- ^)-stretch spanner of Lemma 12.11 on 
all levels below i = H — log(n^). Note that 2^ < diam(S'), and so the longest edge added by 
this construction is of length 0(1/^- ) = 0(diam(S')/en^). So the entire cost of all these 

edges is less than 0(l/e) ■ MST(G)o We then remove from consideration all net-points below level 
i = H — log(n^), set L = i, and focus on the remaining M = H — L = O(logn) levels. Note that 
the aspect ratio of this set is only A = . 

We apply Theorem l4.2l to decompose the remaining graph into all sparse sets in time (ddim log^ n, 

and construct the (1 -|- ^)-stretch spanner of Theorem 13.41 for each sparse subset. Since any point 
of G appears in most logn sets (Theorem I4.2l( iv)) the total time to construct all these 

spanners is (ddim /e)^^'^^™^nlog^ n. Then the stretch between any pair of c-neighbors is (1 -|- j^), 
and it follows from Lemma 12.21 that the total graph stretch is (1 -|- e). 

The weight guarantee follows from Theorem I4.2l fiil in conjunction with Theorem 13.4l fil. □ 

Open problems. Our results suggest several avenues for further research. Gan VFd be reduced to 
^-O(ddim)^ to match what is known for Euclidean space? Also, does the simple greedy hierarchical 
spanner alone yield a light spanner? One can also improve the runtime of our construction. 

^ This observation is due to Shay Solomon, and is similar to one made by Arora in the context of Euclidean TSP 
|Aro98 |. 
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